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We analyze and compare the optimality of approximate and probabilistic universal programmable 
quantum processors. We define several characteristics how to quantify the optimality and we study 
^— ^ ■ in detail performance of three types of programmable quantum processors based on (1) the C-NOT 

gate, (2) the SWAP operation, and (3) the model of the quantum information distributor - the QID 
^vq , processor. We show under which conditions the measurement assisted QID processor is optimal. 

We also investigate optimality of the so-called U-processors and we also compare the optimal ap- 
^ ' proximative implementation of U(l) qubit rotations with the known probabilistic implementation 

^ \ as introduced by Vidal, Masanes and Cirac [ Phys. Rev. Lett. 88, 047905 (2002)]. 
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(N ■ I. PROGRAMMABLE PROCESSORS 

Classical programmable processors are realized as a hardware that perform an operation (called computation) on 
a data register according to instructions (program) encoded in a program register (software). It is one of the central 
issues of computer science whether there exist a universal (classical) processor that performs all possible classical 
transformations of the data register of the size of d bits [1] . Let us, for simplicity, consider only reversible classical 
computation (though the conclusions are valid also for irreversible classical programs). The register composed of d 
classical bits can be reversibly transformed in 2 d \ different ways (permutations), e.g. we consider that on a single bit 
we apply only two programs: the identity and the NOT operation. For a single-bit data register the controlled-NOT 
(CNOT) gate serves as a universal classical processor [the CNOT gate is defined as follows CNOT : j,k — > j, j k 
(j, k — 0, 1)]. That is, if the program register (consisting of a single bit) has the value (is in the state) k = then the 
identity operation is realized on the data bit. Similarly, if k = 1 then the NOT operation is performed on the data 
bit. It is clear that such "control" devices realizing different programs for different bit values of the program register, 
are universal programmable processors also for larger data registers. The size of the program register (in terms of 
the number N of bits) is given by the relation 2 d l = 2 N , because 2 N represents the number of different states of the 
O" 1 register of the size N. As a result we obtain that universal processor for d-bit data register consists of approximately 
N w 2 d (d — 1) bits, i.e. it is exponentially large. 
.,— 1 \ In their seminal paper [2] Nielsen and Chuang showed that the quantum analogue of universal programmable 
processor does not exist, i.e. it is not possible even in principle to design a universal deterministic programmable 
quantum processor. Information about the quantum operation cannot be encoded into the state of arbitrarily large 
program register. In order to realize n unitary transformation of the data register one must use n dimensional 
program register. The same holds also in the classical case, but unlike there, in the quantum case even for a single 
qubit the number of possible programs (unitary transformations) is uncountably infinite. This requires inseparable 
quantum systems as program registers, but such systems are usually excluded from the standard quantum theory. 
The nonexistence of universal programmable quantum device is another example of no-go theorems in quantum 
information processing. 

Let us consider two completely positive maps £ , T given by Kraus operators {Ej}, {-Ffe}, respectively. In addition, 
let us assume that these two operations can be realized with a fixed processor G, i.e. a unitary operation acting on 
the joint data plus program Hilbert space TLd ® "Hp- Denote |He), |Sf) the corresponding pure states of the program 
register, i.e. S[g] = J2j E 3 qE] with E s = (j\G\E E ), where T[q] = J2j F jS Fj with Fj = (j\G\E F ), and {\j)} is 
some fixed basis of TL p . Calculating . E^Fj one derives the following identity 
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EjFj = (E E \E F ) I = cl . (1.1) 



This equation is necessary for simultaneous realization of both operations £, T on the processor G. We remind us 
the ambiguity of Kraus decomposition, i.e. when applying this criterion one must take into account all possible 
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decompositions. A special case is achieved when c = 0, i.e. the encoding program states are orthogonal. In such a 
way any finite number of quantum operations can be realized by encoding them into mutually orthogonal states. Let 
us apply this condition to the case of unitary operations U\, . . . , U n . For them the Kraus decomposition is unique and 
we obtain the set of conditions Ujllk = Cjkl- Obviously Cjj — 1, but in all other cases (j ^ k) it is necessary to set 
Cjk = in order to satisfy Eq. (1.1). It means that encoding n unitarics requires n-dimcnsional program space. By 
construction it can be seen that such dimension of the program Hilbert space is also sufficient. One can simply define 
the processor as G = • Uj®\j){j\. Using this criterion of compatibility (1.1) one can investigate the programmability 
of different families of processes. For instance, in Ref. [3] we have shown that the family of phase damping channels 
can be implemented on a quantum processor, but the family of amplitude damping channels cannot. 

Properties of quantum programmable processors have been studied already by many authors and from different 
perspectives [4-10]. In the present paper we focus our attention on optimality and universality of approximate and 
probabilistic programmable quantum processors. The paper is organized as follows: In the Section II we study 
optimality of programmable quantum processor, in the Section III we analyze in detail optimality of three models of 
programmable quantum processors: the CNOT, the QID and the SWAP processors. The Section IV is devoted to 
approximative programming. We will show in which sense the QID processors are optimal. In the Section V we relax 
the universality condition and we discuss programmable processor that allows us to implement one-parametric group 
of singlc-qubit rotations. We conclude our paper in Section VI with a brief summary of our results. 

II. UNIVERSALITY AND OPTIMALITY 

Even though universal deterministic programmable quantum processors do not exist [2] one can investigate various 
approximations of these processors (this is a general approach when one deals with quantum-mechanical no-go theo- 
rems). One can study scenarios how to achieve the universality by relaxing few of the ideal conditions. In principle, 
there are two options: i) an approximative implementation of quantum programs, or ii) a probabilistic implementation 
of quantum programs. In the first case we allow some imprecision e in the implementation of the desired quantum 
programs, whereas in the second case we relax the condition that the programmability is deterministic by introducing 
a concept of the success probability P SU ccess- Important point is that even though the implementation of the program 
is only probabilistic, the measurement outcomes tell us exactly when the desired operation is performed. It is not 
difficult to see that in both of these cases the universal programmable processors (either approximate or probabilistic) 
do exist [4, 6, 11-13]. However, since the universality is conditioned by some imperfections the question of optimality 
of encoding of quantum operations into states of quantum program registers is of importance [6, 13, 14]. 

For a given processor G one can should study whether it is universal in an approximate, or in a probabilistic sense 
(or both). Different (approximate or probabilistic) universal programmable processors can be compared with the help 
of the approximation parameter e, or the probability success P S ucccss, respectively. There are several ways how to 
characterize the optimality via these parameters. Let remind us the exact definition of these parameters: 

The approximate programmable processors perform a quantum operation T with the precision e(T) = 
minj L>(T, £%), where £(\g] = Tr p G(g £g> £)C?t and D(., .) quantifies the distance between two quantum operations. 
In the probabilistic case we perform a measurement on the program register at the output of the processor. In prin- 
ciple, we can distinguish two cases: either the measurement is fixed (measurement-assisted processor), or the choice 
of an arbitrary von Neumann measurement is a part of the quantum programming. After recording an outcome m 
of the measurement of the program register the data register is transformed into the state g m = T m [g] — j-l m [g], 
where p m = Tr2 m [g] and X m is a linear completely positive, but not necessarily trace-preserving, map. Only in cases 
when p. m p m (g) the transformation T m corresponds to some quantum operation, i.e. it is a completely positive 
trace-preserving linear map. Without the measurement, or better to say without the post-selection, the data register 
is transformed by some quantum operation £ These maps are always performed with a probability equal to unity. 
For each program state one can express the realized operation as a convex combination — Y]j qjTj and the problem 
is whether it is possible to find a measurement M such that the operation Tj was realized with the probability qj . The 
decomposition of £ ^ — q m %n is called realizable if there exists a measurement M of the program register with 
outcomes m such that T m — j-T m . The success probability P SUC ccss of the operation T is defined as the maximum 
of probabilities p over all program states £ with the realizable decomposition £^ = pT + (1 — p)N . The situation 
is simple if one uses the measurement-assisted quantum processor, i.e. the measurement is fixed for all inputs. The 
universality of such device was demonstrated explicitly in Ref. [11], but the optimality is still an open question. It 
is clear that limits of measurement-assisted processors are stronger than limits for probabilistic processors where one 
can vary measurements in order to increase success probabilities. 

There are several approaches how to compare performance of quantum programmable processors. One can use either 
extremal (worst) cases, or average values to evaluate the accuracy /success of approximate/probabilistic processors, 
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i.e. 

Access = J dTP succcss (T), e G = J T dT e (T) 
Success = ™ n Success (T) , e G = max r e(T) 

There is also a freedom in the choice of the function !?(., .) in the definition of the accuracy e(T) of the approximation. 
The larger the success probability the better the probabilistic processor and similarly, the smaller the error the better 

the approximate processor. We will use the notation P^'J^ ess and P su ' ccess for the parameters of the measurement- 
assisted quantum processor specified by the unitary transformation G and the program-register measurement M . The 

following relations hold: max M P S ucccss < ^success and max M Success < ^success- 

Also the universality of quantum processors is usually understood in two different ways: i) either with respect to 
an implementation of all quantum operations, or ii) with respect to a realization of all unitary transformations of the 
program register. In some cases the set of implemented operations can be reduced (restricted) to smaller families of 
quantum operations. It is clear that due to a unitary representation of any completely positive tracepreserving linear 
map both of the meaning of universality are closely related. In order to define optimal quantum processor we usually 
fix the size of the data and program register, d = dimHd and N = dim7i p . The main open problem of quantum 
processor's optimality is to find this functional dependence for success probabilities and approximation accuracy. In 
other words, for a fixed length of data and program register the task is to find the class of processors maximizing the 
success probability and minimizing the approximation parameter. 

Let us take a process fidelity to quantify the distance between two maps, i.e. D(£,T) = 1 — F(£,T) = 1 — /(<&£, 
where $5 = £ ®I[^> + ] ($ + is the projector onto maximally entangled state \ip+) = ^ |j) ® li)) an d .f(Q, a ) = 
(Tr^y/ga^g) 2 is the state fidelity function. From the definition of the quantities P SU cccss (T) , and properties of 
process fidelity it follows that [13] 

Pcrror(T) > e(T) , (2.1) 

where P e rror(^~) = 1 — Psucccss {T) . The probabilistic realization of T means that a given program state £ induces 
the map £^ = P success T + P crlor Af. One can say that program encoded in state £ approximates the transformation T 
with the precision quantified by e = 1 — F(£^, T). Using the concavity of process fidelity one can directly show that 
e < 1 - Psucccss + P C rrorP(A/', T) < 1 - P SU cccss which gives the above inequality. 

Similar relations hold also for the derived average, and worst case quantities. This inequality is not saturated in 
general (for fixed quantum processor), however it might be the case that for the optimal values (optimized over all 
processors) these two numbers coincide. It could be an interesting result if, moreover, the same optimal processors 
are optimal for approximate as well as probabilistic scenario. 

For our purposes it will be useful to have an expression for the process fidelity between an arbitrary channel £ 
(defined as £[g] = J2 r A r gA],) and a unitary transformation U. The state $y = U <g> I[^f+] is pure, i.e. y/Wjj = 
It follows that 

F(U,£) = ($u\£® J[*+]|$ £ ) = ±52V\U*£[\j)(W\k) = ^ITrC/U^. 

j,k r 

III. PROGRAMMABILITY OF UNITARY TRANSFORMATIONS 

In this section we will pay attention to a simpler problem of implementation of all unitary maps. In such case there 
exists a unique Haar measure dU that enables us, in principle, to calculate also the average error and average success 
probability. In what follows we shall analyze three examples of quantum processors. 

A. Controlled NOT 

At the beginning of this paper we have seen that the controlled NOT (CNOT) gate serves as universal classical 
programmable processor implementing the programs on a single classical bit. Its quantum version, Gcnot = I ® 
|+)(+| +er z <g> |— )(— | (with |±) = -^|(|0) ± |1))), is not universal, but still it can serve as a good and simple example of 
a quantum processor realizing approximatively and probabilistically a specific subclass of unitary operations. Using 
this device we are able to realize deterministically arbitrary channel of the form £[g] = pg + (1 — p)cr z g<7 z . Vidal 
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and Cirac [15, 16] showed that using the measurement-assisted CNOT processor one can probabilistically implement 
arbitrary U(l) rotation U v = cnj>(—ipa z ) with the probability P SU ccess(U v ) = 1/2. The fixed measurement is specified 
by the basis |0>, |1). 

In Rcf. [17] we studied different choices of measurement for the CNOT processor. We found that arbitrary measure- 
ment enables us to realize probabilistically a whole set of unitary transformations U v , though in this case probability 
distributions P^ ccesa (U<p) are different for different measurements M. Of course, the transformations I,a z are al- 
ways implemented with the probability equal to. The success probability for a fixed measurement M is given by 
the relation P SU ccess — cos 2 £ cos 2 77 + sin 2 £ sin 2 77, where the angle 77 specifies the choice of the measurement and £ 
represents the states encoding the unitary transformations U v for ip = arccos( -. cos ^ cos v ). Except the 

-^/cos 2 £ cos 2 Ty-t-sin 2 £ sin 2 r\ 

cases when 77 = kir/2 (for k = 0, 1, 2, 3, . . .), each unitary operation {/ is realized. Therefore, either /, or a z is realized 
on measurement-assisted CNOT processor with the smallest probability depending on whether sin 2 77 is larger than 
cos 2 77, or not. In particular P™®^' 71 = minjcos 2 77, sin 2 77}. The average success probabilities for measurement-assisted 
CNOT are given as follows 



Psucc^ = ^ (^J o d £ sin2 £ + cos2 V d£ cos 2£ 



-CNOT.r, 1 ( r W . 2 , , 2 „ „^ 1 

2 ' 

A straightforward calculation shows that in the average sense all measurement-assisted CNOT processors are equiva- 
lent. On the other hand they are not equivalent in the sense of the worst success probabilities. 

It is not easy to calculate these quantities for the case when we are allowed to alternate (optimize) measurements in 
order to increase the success probability. In particular, it is difficult to express analytically the dependence £ = £(</?) 
and consequently to find an analytic expression for the optimal success probability P S uccess(C^)- It is clear that each 
of the transformations can be realized with the probability strictly larger than 1/2, i.e. P S uSSss > 1/2 and also 

-pCNOT 

success I 

Let us assume that we want to perform unitary transformations U v approximatively. What is then the optimal 
program state ^ v approximating the given U v on the CNOT processor? Denoting by £ p [q\ = PQ + (1 — p)a z Q&z an d 
using the derived expression for the process fidelity between an arbitrary channel and a unitary transformation we 
get for the error 

e(U v ) = l-maxF((7 v ,£ f) ) = l-imax{|Tr[/t| 2 ,|TrL/t ( 7 z | 2 }, (3.1) 
P 4 

i.e. it is optimal to encode the operation U v into one of two program states encoding the identity and a z operation, 
i.e. state |±). The encoding is optimal for the state for which the overlap between the desired transformation 
and the program I (or a z ) is larger. In particular, e{U v ) = 1 — maxjeos 2 ip, sin 2 ip}. That is, the characteristics 
of the approximate programmable CNOT processor describing the quality of the realization of the set of unitary 
transformations U v are the following 

ccnot = maxe(l7 v ) = l/2 ) (3.2) 

ccnot = ^Jdpe(U v ) = 1/2-1/tt. (3.3) 

One can directly verify that the bound P e rror > e holds in all its variations. It is saturated only for the worst case 
typical for the optimal measurement-assisted CNOT processor compared with the worst case approximation error, 

, , r ,CNOT,?j=7r/4 , /n 

when we have i-c rr or =£cnot = 1/2. 



B. Quantum information distributor 



Originally, the quantum information distributor (QID) was introduced as a device performing the optimal quantum 
cloning and the optimal quantum NOT operation. In [11] the authors showed that the QID can be used as a 
universal measurement- assisted processor realizing all unitary transformations. The QID is a device with a data 
register representing by a qudit (d dimensional quantum system) and program register represented by two qudits, 
i.e. N = d 2 . It belongs to the family of U processors [3], or equivalently the controlled-U gates. It is defined as 
Gqid = Ek U k ® \0k)(0k\, where k = 0, . . . , d 2 - 1, U k = U ab = £ r e l2 ™ r / d \r - b)(r\ (for qubits sigma operators), 
\6k) = Uk® I\9o), |#o) = \i) ® I n other words, the unitary operators Uk form an orthogonal operator basis 

TrUjUk = dSj k I, and the states \9k) form the basis of two-qudit Hilbert space composed of maximally entangled states. 
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In particular, for qubit Uk = &k arid {|#fc)} is the Bell basis. The program states \9k) encode unitary transformations 
Uk- A general program state £ induces one of the generalized Pauli channels £^[g] — *Ylik'PkUkQU\. withp^ = (9k\£,\Qk)- 
Using the measurement basis \rrij) = \m xy ) = \ — x) <g> J2 r e l2lx vl d \ r — x) one can probabilistically implement 

an arbitrary unitary transformation with the probability P SU ccoss = l/d 2 [12]. The result |moo) = 3? X)fe l^fc) = 
|0) ® -j= ^2 r \r) indicates the successful realization of the program U. In particular, using this measurement basis the 
action of the processor can be written in the following form 

G QID |V> ® |S) = 1 Yl U-jMm}m ® K-) , (3.4) 

3 

where A(S) = 5Zj.(0fc|S)E/fc- If -4(2) is a unitary operation then probabilities of outcomes nij are data independent 
and they read p{rrij) = The described measurement-assisted quantum processor implements unitaries with the 
success probability 

— QID.M = qid,M = , i ,2 ( o r.\ 

success success / ^ ' V 0-, -V 

It is of interest to show whether this is the optimal measurement-assisted probabilistic realization, or not. Note that 
this processor is indeed very specific, because it performs all unitaries with the same probability P S uccess(C^) = l/d 2 . 
A similar conclusions like in the case of CNOT processor can be made. In particular, this measurement-assisted 
QID processor is optimal in the sense of worst case optimality, because using different measurements will result in 
nontrivial success probability distributions over the set of unitary transformations. 

The approximate realization is again similar to the case of CNOT gate. Both of them belong to the family of U 
processors, i.e. they are of the form G — Uj <g> where Uj are arbitrary unitary operations and program states 

{\j}} form an orthonormal basis of H p . These processors arc not very "rich" from the point of view of approximative 
implementation of all unitary transformations. In fact, these processors deterministically perform the random unitary 
channels £p[g\ = J2j PjUjQUj with pj — providing that the program register was initially prepared in the state 

£. As we shall see for the approximative programming only the basis states \j) are useful. In fact 

e(U) = max J F(J7,%) = 1-4 ^Pj\TrU^U \ 2 = 1-4 maxjlTr^^l 2 } . (3.6) 

p d . d J 

3 

For a unitary U the best approximation is achieved by the program state \j) maximizing the overlap \TiU^Uj\ 2 . For 
the QID processor an arbitrary U can be written as a linear combination of the operators Uk, i.e. U = ^2 k ak ^ k anc ^ 
e(U) — l — ^T maxj \a k TTUkUj\ 2 = 1 — max^ | | 2 . Thus, the quality of approximations on the QID processor is given 

by 

eqio = maxe(C7) = 1 - l/d 2 . (3.7) 
This value is achieved for the unitary transformation U = \ Uj . 

C. SWAP 

The quantum SWAP gate acts on two registers of the same size, d = N, and its performance is defined in the 
following way GswAp(g® O^swap = £® B- Taking the SWAP gate as a processor we can implement all contractions 
to a single point specified by the program state, i.e. C^[g] = £ for all input states g. It is clear that such processor does 
not belong to the family of U processors. Indeed it cannot be used to realize any unitary transformation. However, 
unlike for U processors, for each program state the SWAP processor implements different program, which makes it 
exceptional, because there is no redundancy in the state space of the program register. Let remind us that for U 
processors all states having the same diagonal elements in the basis of vectors \ j) encode the same quantum operation. 
The set of contractions is closed under convex combinations and except extremal points of the whole set of quantum 
operations (pure-state contractions) they contain also totally random channel, i.e. a contraction to the total mixture. 

Let us consider an arbitrary measurement M of the program register. The probability of measuring the result |ra) 
is given as 

p(m) = Tr[Q m G SWA p(Q ® £) G swA-pQm] = Tr[£ <g> |m)(m|p|m)(m|] = (m\g\m) , 
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where we used the notation Q m =7® |ra)(m|. That is, for any measurement the resulting probability distribution 
depends on the initial state of the data, i.e. no specific quantum operation is associated with particular results. Or, 
to be more precise, for each outcome the contraction to a fixed point £ is realized. Therefore, measurement-assisted 
SWAP processors are in some sense trivial and do not provide us with any improvement of performance. Therefore , 
it does not make much sense to apply measurements at the output of the SWAP processor. Consequently, the SWAP 

processor is not suitable for performing probabilistically any unitary transformation, i.e. P^^fs — -^success = 0- 

In what follows we will see that from the point of view of approximative implementation of unitary transformations, 
the SWAP processor acts much better. Let us calculate the error e(U) = 1 — ^ max^ ^2 r |TrC/^yl r (S)| 2 , where A r (E) 
are Kraus operators associated with pure program states, i.e. ^4 r (5) = (r|Gg W Ap|2) = (r| ^ . k \k)(j\ <g> \j)(k\E) = 
|S)(r|. In the matrix (basis) representation of Gswap one can use any basis, i.e. even a basis containing the state 
vector |S). Our task is to find e(U) = 1 — y max^ \{ r \U |5)| 2 . Choosing the basis containing the vector |S) 
we obtain that we have to maximize the length of the column of the unitary transformation U'. But we know that 
columns form mutually orthonormal vectors, and therefore in this basis the expression ~^2 r \(r\U' |S)| 2 = 1. This 
means that each pure program state |S) approximates each unitary operation with the same accuracy measured by 
e(U) = 1 — l/d 2 . And consequently 

cswap = £swap = 1 - 1/d 2 . (3.8) 

Thus, the worst case accuracy of the QID and the SWAP processors are the same. Nevertheless, it should be noted 
that for the QID the program space is twice as large as the program space for the SWAP processor. Therefore, we can 
conclude that the SWAP processor is more optimal (suitable) for approximate programming. However, the situation 
is completely different for probabilistic programming. 



IV. OPTIMALITY OF APPROXIMATE PROCESSORS 



Each quantum processor induces a mapping Q from the set of program states into a subset of all quantum operations 
applied on the data register, i.e. Q : £ i— > £^. Let us denote by Tg C L the subset of deterministically implcmcntablc 
quantum programs, where T stands for the set of all possible quantum operations. Since Tg is a linear image of the 
set of program states, it is convex. The question of optimality then can be illustrated in the following way. Denote 
by dT the boundary with respect to some topology. Then the worst case optimality parameter measures the distance 
between the points of the sets Tg and dT. Formally, eg = niaxrer s min^gr D(£,T). 

For the process fidelity we have 

e G = l- min maxF(£,T). (4.1) 

i 61g £61 

Let us consider £' ,T' that optimize the accuracy. Because of the fact that F(^2 k pk£k,1~') > ^2kPkP(£k,T') it 
follows that £' can be always chosen to be the extremal quantum operation of the set T. Let us formulate the main 
problem: given a data register of size d and given the program register of size N. What is the optimal approximate 
processor? This problem is indeed difficult and only partial results arc known. 

Let us pay attention to an optimal realization of all unitary transformations. In this case the approximate processor 
is universal if it performs an arbitrary unitary transformation with nonzero accuracy, i.e. < e(U) < 1. For instance, 
the CNOT processor is not universal in this sense, but the SWAP processor is. The necessary and sufficient condition 
for a processing being universal is the following: The process fidelity has to nonzero for all unitary transformations 
which means that £) — ^2 r \TrU^A r (£)\ 2 > 0. The general processor can be written in the form G = 
^2j k Ajk ® Providing that Ajk form a complete operator basis the process fidelity is different from zero for 

all unitary transformations. In particular, for the SWAP processor the operators Ajk = \k)(j\ undoubtedly form a 
complete basis. 

For the rf-dimensional data register the processor can be universal (implementing all unitaries) only if the operators 
Ajk form an operator basis, i.e. they are independent and the total number of them is d 2 . This necessarily means 
that universal processors for the d dimensional program registers must use at least d dimensional program register. 
It is easy to see that U processors with such program size cannot be universal, because they do not contain sufficient 
amount of independent operators, Ajk = SjkUj. The existence of an approximate universal processor of this size is 
guaranteed by an example of the SWAP processor. The question is whether this processor is optimal, i.e. whether 
eswAP = 1^1 /d 2 attains indeed the minimal value for a universal processor of the program size N = d. D'Arianno and 
Perrinoti [6, 14] found that this is indeed the case for a single qubit (although they considered a different processor). 

In some sense such an error can be achieved trivially. The process fidelity will be nonzero for all quantum operations 
providing that at least for one of the program states £ encodes an operation T such that the state $r has the full rank, 
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i.e. Tr\/ y/Wg^rV^s > for all £. This condition is only sufficient, but not necessary for a processor to be universal. 
Anyway, once we are discussing the problem in general, this condition is sufficient to derive the the lower bound on 
the accuracy of approximate processors. Consider that the processor enables us to perform a contraction into the total 
mixture, i.e. let us consider the map A[g] = \l ■ Then for any £ the process fidelity reads F(A,£) = 37 [TrV*J] 2 , 
because $_4 = -^1. The minimum is achieved for pure states <E>£ corresponding to unitary transformations. That is, 
all other programs are approximated by this operation better than unitaries and therefore for such universal processor 
e<3 > 1 — 1/d 2 . The SWAP as well as the QID processors can realize the transformation A and they both saturate 
this bound. 

Let us restrict the original problem and ask the following question: What is the optimal U processor? To be 
able to discuss optimality processors they have to be universal. It means that this question does not make any 
sense for program registers of the size N < d 2 . The QID is an example of an approximate processor in the case 
when N — d 2 . For the U processors the situation is simpler because the process fidelity is given by the formula 
F(U,£) = -^J2 r |Trl^A r | 2 . Let us fix the dimension of the program space to be N = d 2 , i.e. we work with d 2 
unitaries Uj. The first question is which of the U processors is optimal in implementing all unitary transformations. 
We know that the best approximation is achieved for the basis states \j) and e(U) = 1— -h m&Xj{\TrU'Uj\ 2 }. That is, 
the worst case is represented by a unitary U x having the smallest maximal overlap with "elementary" programs Uj . Our 

aim is to show that this worst case is optimized for mutually orthogonal operators (jj, i.e. TrcrJ Ok = dSjk- Processors 
having as elementary programs mutually orthogonal unitaries are all equivalent. In particular, the approximation of 
the operator U x = g ° j maximize the error to e(U x ) = 1 — 1/d 2 . Each unitary transformation can be expressed in 
a suitable orthogonal operator basis in such form, i.e. each unitary operation can represent the worst case for some 
processor. 

Consider that we know the transformation U x for a processor with nonorthogonal unitaries. Let us define a processor 
using orthogonal operators Uj such that U x = ^ J2j a i ■ Let us express also the elementary operators Uj in this new 
operator basis Uj = J2 a u ja&a- Calculating the overlap we obtain the following bound 

\TrU x U k \ 2 = \^u kj \ 2 < Yl K-l 2 < 1 . ( 4 - 2 ) 
j j 

where we used the identities Trover;, = d8 a b and Tr/T^l^ = dJ2j \ u kj\ 2 = d. Consequently, the error is minimized for 
I Tr a*; 1 2 = 1 and therefore we can conclude that the U processors with orthogonal elementary programs are indeed 
optimal, i.e. the QID processor is an optimal approximate U-processor for TV = d 2 dimensional program register. 

This optimality result for approximate processors can be used directly for probabilistic processors by using the 
inequality P rror > e, or equivalently Psucccss < 1 — e. The problem is now the following: Which of the U processors 
is optimal in a probabilistic realization of unitary transformations? Based on the answer to the similar question 
for approximative programming we can say that the optimal success probability is less or equal to 1 — e = 1/d 2 , 
i.e. P S ucccss < 1/d 2 . However, we know that the measurement-assisted QID processor saturates this bound, i.e. 
^successf — 1/cZ 2 , and therefore the QID is an example of the optimal probabilistic U processor implementing all 
unitary transformations. 

V. OPTIMALITY FOR U(l) ROTATIONS 

In this section we will relax the universality condition and we will study an implementation of a one-parametric 
group of unitary transformations. An example of the universal processor performing this task is the CNOT processor. 
Our aim will be to specify the dependence of the approximation error and the success probability on the size N of 
the program register. 

Let us start with the approximative realization of U v — exp(i(pA) on U processors, i.e. we use processors of the form 
G = Uj ® \j}(j\, where j = 1, . . . , N. Since the approximation error is specified by overlaps between U v and Uj, 

it is reasonable to assume that Uj are from the linear span of unitaries U v . Otherwise, we would obtain TrU^Uj = 
which is not interesting for approximative realization of U v , i.e. the associated states \j) do not approximate any 
transformation U v . 

In particular, Uj are from the group U v . The overlap between two unitaries U V ,U V ' is given by the expression 
iTr/TjLVl 2 = |Trexp[i((/ - Lp)A]\ 2 = \ J2 n e^'-^" | 2 , where a n are eigenvalues of A. The closer the angles <p',<p 
are, the larger is the overlap. In what follows we will assume that the data register is two-dimensional, i.e. A = a ■ a, 
U v = cosi^/ + isinif(a ■ a) and |Tr/7|[/^/| 2 = 4cos 2 ((^' — <p). 

Let us consider a processor with dimensional program register performing TV unitaries U v . The question is, 
what are the best choices of angles (fj in order to minimize the approximation error of the implementation of the 
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whole group U v {<p G [0, 2ir]). The program state \j) best approximates those angles ip that lie in the vicinity of 
the angle ipj (the difference ipj — ip is smaller than the differences ip^ — tp for k j). Because of the identity 
cos (<Pj — f + 7r) = cos 2 ((/9j — ip) the state \j) approximates with the same accuracy the unitaries U v and U v +k- That 
is, it is enough to consider only implementation of unitaries specified by angles within the interval ip £ [0,7r]. The 
whole problem can be illustrated in the following way. The angles ipj divide the half-circle (angles from to n) into N 
regions. The state \j) approximates optimally all the angles from the interval [tfj — \{<Pj — <Pj-i), <Pj + 5(^+1 — <Pj)] 
(for simplicity we assume that tpi = and ipj < <Pj+i). The best choice of tpj is such that the half-circle is divided 
into equally large regions, i.e. the angles ipj are separated by the same angle tpj+i — <Pj = n/N. In such case the 
approximation error reads 

e G = 1 - -U cos 2 (tt/ (27V)) = 1 - cos 2 (7r/(2iV)) (5.1) 

because for arbitrary <p the smallest difference is minj{pj — ip} < n/(2N), i.e. the overlap is maximal and the error 
£g is optimal. The average error of implementation of the operations U v equals to 

Based on this result we can bound the optimal probabilistic realization of qubit operations U v using the N dimen- 
sional program register and the U processor 

Success < COs2 (^)' ( 5 - 3 ) 
q 1 N TT 

^success < 2 <yl + ~ S[ll N^ ^'^ 

Now the question is whether there exists a measurement-assisted U processors saturating these bounds. In what 
follows we will present an example of the U processor that saturates this bound in a limit sense, i.e. for large program 
registers. 

In Rcfs. [15, 16] the authors proposed a way how to increase the success probability to arbitrarily close to unity 
by using large program registers with the dimensionality N = dim7i p = 2", where n is the number of qubits in the 
program register. They found a network of the controlled U form composed of n sequentially applied fc-Toffoli gates 
for k = 1, 2 ... n. The fc-Toffoli (Tk) implements a controlled NOT operation with k control qubits and a single target 
qubit. The whole processor acts onn + 1 qubits, where the first qubit represents the data and all other qubits form 
the program register. The fc-Toffoli gate uses qubits 1, . ..k as control qubits (this includes the program qubit as 
the qubit number 1) and the (fc + l)th qubit is the target qubit. The whole processor is described by the operator 
G = T„(T„_i ® I) . . . (T 2 ® 7®(™- 2 ))(Ti ® Z®^" 1 )). Vidal et al. showed [15, 16] that it is possible to choose program 
states and a measurement such that the success probability scales as 

PsucceM = 1 - (1/N) (5.5) 

for an arbitrary operation U v = e lvrTz . Thus, we find that P^'Jtf css — -f^lLss = 1 ~ (V-^Oi which is consistent with the 
derived bound, but scales differently, because P^ uccess < 1 — sin 2 (7r/27V) — > 1 — (n/(2N)) 2 for large N. This means 
that the saturation of the bound is still an open question. Moreover, let remind us that the derived formula for the 
measurement-assisted U processor holds only for dimensions N = 2", where n is the size of the program register in 
the number of qubits. 



VI. CONCLUSION 



One of the goals of the research in the field of quantum computing is a construction of programmable quantum 
processor. It follows from the work of Nielsen and Chuang that perfect programmable processor does not exist. 
In this paper we have analyzed two different scenarios relaxing the condition of universality: approximate imple- 
mentation of quantum programs (approximate processors) and probabilistic implementation of quantum programs 
(probabilistic processors). We have discussed the problems of optimality of universal approximate processors and 
universal (measurement-assisted) probabilistic processors. Using examples of the CNOT, the QID and the SWAP 
processors we have shown explicitly the validity of the general relation between the approximation error e& and the 
success probability P B „ ccess , respectively. In particular, we have seen that the inequality P crr0 r > e is not saturated in 
general case. For instance, the SWAP processor is completely useless in probabilistic programming P s r ^ p = 1, but 
approximatively it performs optimally. 
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We have studied in detail optimality of restricted class of processors - the so called U processors. We have shown 
that the QID processor is optimal in implementing unitary transformations on a qudit with N = d 2 dimensional 
program register. Under the same settings we have analyzed the restricted universality of an implementation of only 
one parametric set of unitary transformations. Unlike the QID processor we have found that the presented example 
of the U processor does not saturate the bound derived from approximative implementation. That is, the question of 
optimality for this processor is an open question. 

The optimality questions of universal processors (either approximate, or probabilistic) represent a difficult open 
problem of quantum information science [18]. In this paper, we have described current state of the art and we have 
presented known results. The question of efficient programmability of quantum computers makes these problems of 
optimality very attractive and they deserve further investigation. 
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